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Representations  based  on  Zero-crossings 
in  Scale-Space 


Robert  A.  Hummel 


Abstract 


Using  the  Heat  Equation  to  formulate  the  notion  of  scale-space  filtering,  we  show 
that  the  evolution  property  of  level-crossings  in  scale-space  is  equivalent  to  the  max- 
imum principle.  We  briefly  discuss  filtering  over  bounded  domains.  We  then  con- 
sider the  completeness  of  the  representation  of  data  by  zero-crossings,  and  observe 
that  for  polynomial  data,  the  issue  is  solved  by  standard  results  in  algebraic 
geometry.  For  more  general  data,  we  argue  that  gradient  information  along  the 
zero-crossings  is  needed,  and  that  although  such  information  more  than  suffices,  the 
representation  is  still  not  stable.  V/e  give  a  simple  linear  procedure  for  reconstruc- 
tion of  data  from  zero-crossings  and  gradient  data  along  zero-crossings  in  both  con- 
tinuous and  discrete  scale-space  domains. 


1.   Scale-space  and  zero-crossings 

The  use  of  multiresolution  representations  is  an  important  idea  for  the  analysis 
of  signal  and  image  data.  Many  data  structures  have  been  studied,  including  Gaus- 
sian pyramids,  difference-of-Gaussian  pyramids,  Laplacian  pyramids,  and  "scale- 
space"  formulations  [1,2,3].  The  latter  formulation,  to  be  described  briefly  below, 
can  be  used  as  a  continuous  model  of  the  other  formulations.  We  will  discuss  the 
representation  of  data  by  zero-crossings  in  scale-space,  and  consider  the  stability  of 
reconstruction  methods. 

The  natural  framework  for  the  analysis  of  scale-space  formulations  of  mul- 
tiresolution representation  is  in  terms  of  the  heat  equation  [4,5].  Specifically,  let 
fix)  be  a  bounded  function  defined  for  x^lR".  (Arbitrary  dimensions  can  be  han- 
dled with  little  additional  fuss  over  the  case  n  =  1;  we  will  later  comment  on  the 
case  of  a  bounded  domain  D  C  IR"  ).  We  define  u(x,t)  to  be  a  bounded  solution  to 
the  heat  equation: 

du         . 

—  =  Am,  (Heat  Equation) 

u(x,0)  =  fix). 

The  solution  is  given  by  convolution  against  the  fundamental  solution  to  the  Heat 
Equation,  which  for  the  domain  IR"  is  given  by 

uix,t)  =    SKix-y,t)fiy)dy, 
IR" 

where 
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Kix.t)  =  (4'iTr)-'"'2e-kP/4f 

We  see  that  u(x,t)  is  obtained  by  blurring  f{x)  by  increasingly  diffuse  Gaussians, 
parameterized  by  r>0,  with  standard  deviations  ct  satisfying  2ct^  =  4r.  In  computer 
vision,  scale-space  sometimes  refers  to  the  (x,cr)  variables  that  can  be  used  to 
reparameterize  the  domain  of  u.  We  retain  the  {x,t)  parameterization  to  keep  the 
linear  Heat  Equation  relation  for  the  function  u. 

Convolution  by  Gaussians  is  considered  special  for  many  reasons  [6,4,7].  We 
see  from  the  above  analvsis  a  relationship  betv/een  Gaussian  convolution,  the  Heat 
Equation,  and  the  Laphcian  operator.  „Gf  course,  Gaussian  convolution  enjoys 
other  properties;  for  example,  ihe  central  Jjmit. theorem  implies  that  Gaussian  con- 
volution is  easy  to  implement  by  an  iterative  procedure.  However,  we  also  see  the 
extent  of  the  similarity  of  difference  of  Gaussians  and  the  Laplacian  of  the  Gaus- 
sian; namely,  since  Kix,t)  is  itself  a  solution  to  ^qjlieal;  Equation, 

(AA)  Kx,<r^"f(Jc;if)  ^«^(xv?M't)  -  Kix,t))/r. 

That  is,  the  difference  of  Gaussians  is  a  good  approximation  to  AAT  as  the  separation 

between  the  spread  of  the  t\vo  Gaussians  approaches  zero  (and  the  difference  is 

scaled).  "    *'     '  '  ' 

'  jfT.'O':  '-.rr.'- 
Filtering  by  the  Laplacian  of  a  Gaussian  can  b*^  written  in  three  ways: 

If  we  denote  the  result  by  v(xiif),-we^e  l5tat 

(1)  v{x,t)  is  the /(x)  data  filtered  by  the  Laplacian  of  a  Guassian. 

(2)  v{x,t)  is  the  solution  to  the  Heat  Equation  with  initial  data  A/. 

(3)  v{x,t)  is  A«(x,r),  where  u  is  the  solution  to  the  Heat  Equation  with  initial  data 
fix). 

The  zero  set  of  v{x,t)  is  the  point  set  in  {x,t)  where  v  =  0.  The  set  might  be 
empty  (for  instance,  if  /  is  subharmonic  or  superharmonic;  [8])  everything  (if  /  is 
harmonic;  [9]  ),  or  a  proper  subset  of  {x,t)  space.  In  the  latter  case,  zeros  can  be 
isolated  points,  lines,  and  surfaces  (but  never  regions).  We  distinguish  components 
of  the  zero  set  which  form  manifolds  of  codimension  one: 

Definition:  The  zero-crossings  of  v(,x,t)  refers  to  the  point  set 

d{{x,t)\v{x,t)  <Q}    n    d{{x,t)\v{x,t)>Q}.    ■ 

Zero-crossings  have  been  suggested  for  segmentation  of  imagery  by  edge  detection 
[10],  and  for  stereo  matching  and  motion  correspondence  between  pairs  of  images 
(e.g.,  [11]).  It  has  also  been  suggested  [7]  that  the  zero-crossings  are  a  nearly  com- 
plete representation  of  A/.  Finally,  Wilkin  [3]  observes  that  zero-crossings  in 
scale-space  evolve  as  t  increases,  and  are  never  created  at  some  nonzero  t  .  This 
property,  discussed  in  [7]  and  in  [6],  ensures  that  zero-crossing  surfaces  are  nested, 
one  within  another,  enclosing  regions  containing  the  face  {t  =  0},  or  forming  a  sheet 
meeting  the  face  [t  =  0}  and  extending  to  r  =  0°.  The  property  can  be  given  a  pre- 
cise statement: 

Evolution  property  of  zero-crossings:  Let  C   be  a  connected  component  of  the  set  of 
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zero-crossings  in  the  domain  {(,x,t)\x^]R",  Ti^t^Tz},  where  O^Ti  <  Tj-    Then 

c  n  {(x.ok  =  ri}#0.  ■ 

In  the  remainder  of  this  paper,  we  wish  to  make  two  main  points.  First,  we 
establish  the  equivalence  of  the  "Evolution  property  for  zero-crossings"  and  the 
classical  maximum  principle  for  parabolic  partial  differential  equations,  thereby 
allowing  us  to  consider  bounded  domains  and  nonstationary  convolutions.  Second, 
we  consider  reconstructibility  from  zero-crossings.  The  principle  result  is  that  if 
knowledge  of  the  location  of  the  zero-crossings  is  sujxplemented  with  gradient  infor- 
mation of  V  at  only  those  points  in-ihe  zer^-vtdssihgs',  then  there  is  a  simple  scheme 
for  reconstructing  some  of'ih^'^hxA.['~^Mt\^ii  even'th'e-n  numerical  accuracy  of  the 
reconstruction  is  unstable.^^'-^^"'^  ^-/i!.^:t:i  :-f  vo  ins.r.air  -n  c: 

jp.fi  an/;i82i;iO  lo  aonsiSiltb  "io  \tn£.[irr.rc 
2.  The  Maximum  PfMiipIid  o:  norjuioz  £'-.i3^'.:f  ?!  (i.x)^  aonip 

The  classical  maxjiroim.^  pj:jf>ciplfi=  ^91^  )tS^  sali;t^iqn  to  the  parabolic  equation 
duldt  =  A«  states  (see,  e.g.',  [12,16,14]):'      •'^ 

Maximum  Principle:  Let^P  Q  Ig."  be  open  and  bounded.  Suppose  u  is  a  solution  in 
T  =  {(x,t)  \  x^D  ,  0<t<T}  of  class  C^  which  is  continuous  in  the  closure  T.  Then 
u  assumes  its  maximum  at  some  point  {x,t)  for  which  either  x^^P  or  r  =  0.  ■ 

Next,  suppose  that  scale-space  construction  is  denoted  by  the  operator  v  =  Sg, 
which  is  to  say  that  the  scaled  fahctid'nV^,^)  is-obtained  from  the  initial  data  g(x). 
In  the  previous  section,  we  d^ine^..^^  -to  be  v(x,r)  =  K(-,t)  *  g,  where 
gix)  =  ^f{x),  but  we  can  imagine  mpre  general  operators.  In  any  case,  it  is  logical 
to  make  certain  assumptions  ab&ull'S,  although  all  we  will  require  is  that 

(1)  If  ^(a:)  is  continuous,  then  i'^' Is 'Continuous. 

(2)  S{  — g)  =  —Sg  for  all  g,  and  if  v  =  S^  and  g(x)  =  v(x,7),  then  v  =  Sg,  where 
v(x,t)  =  v(x,t  +  T). 

(3)  If  ^(x)-0  as  \x  |-oo  and  v  =  5^,  then  for  each  t,  v{x,t)-4)  as  |x  |-oo. 

We  note  that  if  S  is  defined  by  convolution  with  Gaussians  as  in  Section  1,  then 
the  maximum  principle  holds  for  v  =  5^  as  long  as  g  is  continuous;  further,  condi- 
tions (1)  -  (3)  hold. 

Our  first  result  is: 
Proposition:  The  following  are  equivalent: 

(i)  The  maximum  principle  holds  for  solutions  v  =  Sg  using  continuous  initial  data 
g(j:)  satisfying  ^(x)-0  as  |x|-oo. 

(ii)  The  evolution  property  holds  for  level-crossings  of  solutions  v  =  Sg  using  a 
scale-space  operator  S  satisfying  (1)  -  (3)  above  and  continuous  gix)  satisfying 
gix)^  as  \x  1-00. 

Proof:  We  first  show  that  the  maximum  principle  implies  the  evolution  property. 
For  if  the  evolution  property  fails  for  a  level-crossing  /  for  some  v(x,t),  then  by 
suitably  transforming  v,  g,  and  /,  we  can  assume  that  /^O,  and  that  there  is  a  solu- 
tion V  =  Sg,  (with  ^-O  as  |x|-co)  with  a  component  C  of 
{(x,t)  I  0<r<7,  v(x,r)>/}  disjoint  from  the  plane  t  =  0.  Let  (xo,to)  be  a  relative 
maximum   in   C  .    Then   since   C    is  open,  there  is  a  bounded   cylinder   in   C    with 
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(jco.'o)  in  the  interior,  which  is  in  violation  of  the  maximum  principle. 

Conversely,  if  the  maximum  principle  fails,  then  for  some  cylinder  Px[0,7] 
and  some  v(;c,r)  given  by  v  =  5^  with  gix)^  as  \x  |-oo,  the  maximum  of  v  occurs 
either  in  the  interior  or  the  top  of  the  cylinder.  Either  way,  there  is  a  value  /  less 
than  the  maximum  but  greater  than  the  values  on  the  bottom  and  sides  of  the 
cylinder.  Thus  there  is  a  component  of  the  level-/  crossing  within  the  wedge  O^f  ^7 
which  lies  entirely  within  the  interior  of  the  cylinder,  and  thus  does  not  meet 
{t  =  0}.    So  the  evolution  jprope..^  Tov  level-crossings  is  violated.  ■ 

;<-,  "ic  i!>i::oi)<>"i3]';.'    •.. 

We  illustrate  the  utility  of  tlie  proposition  with  three  observations.  First,  given 
the  equivalence  with  the  maximum  principle,  any  proof  of  the  evolution  property 
for  zero-cro5ding3'-tK^t'3<^§'"ribt'4eftHer''Q§^'1flife'^tfpiximum  principle  or  essentially 
redo  the  proor  of  ihc -m&3efiiiu*S^  p'rinci^i^^i^  hlgfily  suspect.  Since  the  maximum 
principle  is  slightly  di'ilicktiViTcspfeclalTy^j^n'^fiie'Sb^ft^e^^  strong  regularity  assump- 
tions, the  former  course  secmsmoTii  kppropru.te.^   ■'■p 

Second,  using  the  version  of  the  gradient  I^opf  maximum  principle  for  the  Heat 
Equation  [13],  it  is  not  hard  co  show  that  knowledge  of  the  zero-crossings  together 
with  gradient  information  along  v/iih^the  zero-crossings  (or  even  just  one  zero- 
crossing  contour)  is  sufficient  to  determine" v(x,i)  uniquely  [15].  In  section  3,  we 
give  a  more  constructive  diacussion  of  fKis  point;  but  it  is  interesting  that  the  max- 
imum principle  establishes  this  uiliqufehx;s&'}     - 

Finally,  we  see  that  any  scaling  method  obeying  the  maximum  principle  will 
yield  the  evolution  property  for  zero-crossings.  Under  fairly  severe  restrictions, 
this  leads  one  to  Gaussian  convolution  16],  but  more  general  scaling  methods  are 
possible.  For  example,  blurring  tiy  a  pa/abolic  operator  of  the  form  du/dt  =  Lu, 
where  L  is  a  uniformly  elliptic  linear  second  order  differential  operator  with  non- 
constant  coefficients  will  certainly  still  give  a  maximum  principle.  In  fact,  L  can  be 
nonlinear  [16].  Moreover,  suppose  we  replace  IR"  with  a  bounded  domain 
P  C  IR",  and  insist  on  data  f(x)  with  compact  support  in  the  interor  of  P  .  We  may 
then  define  v  =  Sg,  where  g  =  A/,  by  solving 

dv/dt  =  Av    in  P  x(0,oo), 
v(x,0)  =  g{x)    for  x^P 

vixj)  =  0    for  x^dP  . 

This  scaling  is  not  given  by  convolution  against  a  Gaussian,  since  the  domain  is 
bounded,  but  nonetheless  obeys  a  maximum  principle,  and  gives  the  same  evolution 
property. 

3.   Completeness 

We  return  to  a  consideration  of  zero-crossings  of  data  v(x.t)  obtained  by  filter- 
ing initial  data  f  (x)  defined  for  ,r  ^IR"  by  the  Laplacian  of  Gaussians,  AA'(-,f).  The 
question  we  wish  to  address  is:  to  what  extent  do  the  zero-crossings  represent  /(.r)? 
Clearly,  f(x)  can  at  best  be  reconstructed  to  within  an  arbitrary  additive  harmonic 
function  and  a  scalar  multiple.  However,  if  wc  assume  that  f(x)^  as  |.t  |-oo.  then 
only  the  multiplicative  constant  is  of  concern. 
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Yuille  and  Poggio  [17]  make  the  observation  that  if  g(x)  and  hence  f(x)  is 
polynomial  in  x,  and  if  n  =  l,  then  reconstruction  from  zero-crossings  is  theoretically 
possible.   They  also  refer  to  the  validity  of  the  observation  for  larger  n. 

We  note,  however,  that  when  g(x)  is  a  polynomial  in  at^IR"  for  any  n,  then 
v(x,t)  is  a  polynomial  in  (x,r)  ^IR""^^.  Accordingly,  the  zero-crossings  are  part  of 
the  analytic  varieties  of  the  polynomial  v  as  studied  in  algebraic  geometry.  It  is  well 
known  that  the  varieties  in  C  determine  the  complex  polynomial  defined  on  n  com- 
plex variables.  It  is  not  as  commonly  ppj4^-iibi\^  npnetheless  true,  that  an  n- 
dimensional  subportion  of  the  intersection  of  the  analytic  variety  with  R""^^  also 
determines  the  polynomial  ,[18].  Thus  the  case  of  polynomial  data  can  be  settled 
with  algebraic  geometry.  ,   .     ^     ^ 

However,  since  th&_^^Tn^^Qn-^Q^i^^p^y.nQsmh^yiMoY^f^^^^^^  is  essentially 
an  analytic  continuation  r^JfiiiVi  m^^Jf&nft^j^^jR^'^^S^^i'.P^'o^  is  unlikely.  That  is, 
small  errors  in  measurM5i,e5$ogf^t|]efjiEei;9c(g'a5Smg§;  coul4  lead  to  arbitrarily  large 
errors  in  the  determination  o4jg(i>r)&7  Icfi?  ^^If^rently,  there  can  be  widely  different 
initial  data  leading  to  nearly  identical  zero-crossing  data. 

r"  • ; ,     1 '       J  t:  1  .J 

Worse,  settling  the  case  fpr, .polynp'mifil  data  says  little  about  the  case  of  con- 
tinuous initial  data.  Althoi^g^.^he  Stone7^y^ierstrass  theorem  says  that  a  continuous 
function  can  be  uniformly  .approximated,  .bj^j^  polynomial  on  a  compact  set,  the 
zero-crossings  depend  on  the  init;%l  4aitat$i9^^fe»-*'^<^  ^^^  dependence  can't  be  local- 
ized. Further,  the  lack  of  stability  r^^^^s^^-lihi^jtbie, approximation  is  irrelevant.  The 
situation  is  similar  to  the  fact  that  a  polynomial  of  a  single  variable  with  all  real 
roots  is  determined  by  its  zero's  .^^iitth'af  given  all  the  zeros  of  a  continuous  func- 
tion, one  knows  nothing  more  than  the'' zeros.  ' 

In  fact,  there  are  known  exarnples  of  pairs  of  functions  fi(x)  and  f2(x)  such 
that  the  corresponding  vi{x,t)  and  V2ix,t)  have  identical  zero-crossings  at  all  levels 
of  resolution.  John  Daugman  supplies  the  example  (for  two  space  dimensions)  of 
fi(xi,X2)  =  sinxi,  and/2(xi,A:2)  =  (sinxi)(2  +  cosjr2). 

However,  if  the  zero-crossing  data  is  supplemented  with  knowledge  of  the  gra- 
dient data  at  the  zero-crossings,  then  reconstruction  of  at  least  some  of  the  data 
g(x)  is  theoretically  possible  by  a  quite  easy  procedure,  given  below.  Details  of 
these  ideas  were  reported  earlier  in  an  unpublished  work  [15].  The  use  of  gradient 
data  for  the  representation  also  appears  in  [4],  but  the  gradient  data  there  is  not  lim- 
ited to  the  zero-crossings.  The  use  of  gradient  data  along  zero-crossings  is  dis- 
cussed in  [19].  Many  researchers  have  noted  from  a  casual  observation  of  zero- 
crossings  of  image  data  that  zero-crossings  with  large  gradient  magnitudes  are  of 
greater  significance  than  those  with  low  gradient  magnitudes. 

3.1.   Continuous  Case 

Specifically,  let  O  be  a  hounded  connected  component  of 
{(x,t)\  r>0,  v(x,t)=j^O},  and  denote  by  P  the  set  (x  ^IR"  |(x,  0)  ^H},  and  by  T  the 
zero-crossing  afln{/>0}.  Let  t  be  a  value  such  that  T>sup{r  [(x.r)  €0}.  Next,  we 
set  g(x)  -  gix)  for  x6P,  and  g(x)  -  0  elsewhere.  Finally,  let  h(x)  be  the  g(x) 
data  blurred  to  the  level  t: 
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b{y)  =   SKiy-x,'r)g(x)dx. 

Using  Green's  theorem,  it  is  easy  to  show: 

Proposition: 

b(y)  =  SK(y-x,7-t)Vv(x,i)nd<T, 
r 

where  n  is  a  surface  normai  to  f^at  (x,r),  and  da  is  surface  area  measure.  ■ 

Thus  given  the  zero-crossi:  g  Fyand  f^vXixvf),:fi-r  {x,t)^r,  then  the  blurred  data 
b(x)  can  be  constructed  by  a  simple  linear  process.]  The  original  data  g(x)  can  be 
reconstructed  by  deblurring  the  b(x)  data.[iO].  Peblurring  is,  of  course,  a  classic 
unstable  proc.s:  TT^^e- skjiatioo. is^Tjot'-hojjeTessi  how^-iver,  since  g(x)  has  known 
compact  support,  v/hich  might  be  used  to  advantage,  and  also  since  errors  that  occur 
are  predominantly  in  high  frequency  components,  which  might  not  be  as  essential  to 
visual  interpretability.  -   _Lr.iJ-idZiZi_: — 

The  lesson  of  this  section,  ultimately,  is  that  even  for  bounded  zero-crossings 
supplemented  with  gradient  c'ata  along  the  zero-crossing,  reconstruction  is  still 
unstable.  We  defer  a  remark  on  relaxing  the  constraint  that  the  zero-crossing  be 
bounded  until  the  next  subsection,  where  we  "consider  a  discrete  version  of  the  result 
of  this  section. 

3.2.  Discrete  Data  ^  " 

For  simplicity,  we  treat  the  case  of  one  unbounded  space  dimension,  although 
the  results  extend  easily.   We  a/.e  given  d'.ta/i,  i-     •     ,  —  1,0,1,  •  •  ■  ,  and  define 

gi  =  ~^ji-\--:^ji^  4-''  +  i- 
We  define  the  filtered  data  v,-^;^  recursively: 

v,,o  =  gi. 

We  also  define  the  blurring  kernel 

Both  V  and  K  satisfy  a  discrete  version  of  the  Heat  Equation,  namely 

4  2  4 

It  is  not  hard  to  prove  a  discrete  analogue  of  the  evolution  property  for  zero- 
crossings.  The  key,  as  one  might  suspect  from  Section  1 ,  is  a  discrete  version  of  the 
maximum  principle,  which  is  easy  to  establish. 

Let  fl  be  a  bounded  4-connected  collection  of  pixels  {i,k)  with  a  nonempty  set 
D  =  {i|(i,  0)^n}.  Let  T  be  an  upper  bound  T>ma.x{k\(i,k)(.D,},  and  define  h,  to 
be  the  data  gi,  iiD  ,  blurred  to  level  T: 
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Finally,  let 


a(±i,o)fi  =  {u,k)^n\(i±i,k)ia}, 

3(0.1)"  =  {ii,k)^^Hi.k  +  l)i^}, 

3(0,-1)0  =  {{i,k)(:fl\  ^>fi^iLk'-l)i^}- 

Then  simple  but  messy  algafiriiailjovSfe  hsitdlshowcO.:-:  k, 

Proposition:  T    .zzaooiq  -rBsni!  3iqrm?.  ^  vc  jLJ-j.-'an 


46;  =        2  S^Off     ?OTbmTTf--^{^J+«ft9/t«tW-?l-5/n"-j,?^--l] 

e=-l,l    (»,*)€3(W))fti6GJf7S^n)B  01  bsp.u  CG  iri^irn  djiuv.- 

-':.iriw  .ElpsnoqiTOS  v^^rio;;^.!!  ij;  ■' 

^i+e-j.r-t-i+^i-j.r-t'-i ,  , 

^ r [v,  +  e,k-Vi,;fc] 

.  '  j  .""-rn  fir  //IsJfiTii.  !■        ■' 

(<-,*) €3(0,1)0  .1-       -J  IV  3Tcn/  :. 

S        4v,-,f/i:f_j.r-t. 

('.*)€  3(0.-1)" 

To  reconstruct  data  by  the  ab^y^  ^nation,  choose  a  connected  component  of 
{(i,k)\v(i,k)>0},  (or  respectively  <0  ).  If  the  component  extends  to  infinity  in 
either  coordinate,  truncate  the  domain  to  become  a  convenient  bounded  collection  of 
pixels,  and  denote  the  result  by  fi.  We  store  the  sets  3(  +  i,o)0.  3(0,±i)O.  and  D  as 
defined  earlier.  For  pixels  (i.k)  in  3(i,o)"  (respectively  d(-i^o)tl  ),  we  store  the 
information  v,-,jt  and  v,  +  i,jt,  (respectively  v,-,jt  and  v,_i,jt).  For  pixels  ii,k)  in  3(o,i)i^. 
we  store  the  data  v,_jt  +  i  and  for  3(o,-i)n  pixels  we  store  v,_jt.  Using  the  above 
equation,  we  choose  a  T  and  reconstruct  the  blurred  data  bj.  To  reconstruct  the 
data  gi  for  i^D  ,  it  suffices  to  deblur  the  fe,-  data  by  solving  for  gi  in  the  linear  equa- 
tions defining  &,-.  In  fact,  the  system  is  overdetermined,  although  still  poorly  condi- 
tioned, especially  if  \P  |  or  T  is  large. 

In  order  to  make  the  computations  feasible,  it  is  necessary  to  modify  the  for- 
mulas for  a  bounded  spatial  domain.  We  in  fact  solved  a  bounded  domain  problem, 
with  -A^</<A^,  setting  v,-,jt  =  0  for  i  =  ±N.  The  blurring  kernel  K  is  changed  by 
this  modification,  but  the  proposition  carries  over  with  little  change. 

In  Figure  1,  we  show  a  1-D  signal  gi,  and  the  zero  crossing  separating  positive 
and  negative  regions  of  the  associated  v,_jt.  Applying  the  above  procedure  to  the 
central  positive  component  Cl,  we  obtain  the  reconstructed  hi  data  shown  in  Figure 
2.  The  true  i>,-  data  is  identical  to  essentially  machine  precision.  Finally,  using  the 
method  of  pseudoinverses  to  deblur  the  data  shown  in  Figure  2,  we  obtain  gi  for  /  in 
the  middle  range,  as  shown  in  Figure  3.  This  is  to  be  compared  with  the  true  initial 
data  in  Figure  la.  The  poor  correspondence  is  due  to  the  fact  that  the  deblurring 
problem  is  poorly  conditioned;    better  deblurring  results  are  obtained  if  the  amount 


Page  7 


Representation!  based  on  Zero-crossings  In  Scale  Space 


of  deblurring  is  very  small.  However,  this  requires  that  the  top  of  zero-crossing 
contour  enclosing  the  data  occurs  after  not  too  many  blurring  steps.  The  lesson 
learned  here  is  that  although  reconstruction  is  in  theory  possible,  practical  recon- 
struction may  be  impossible. 
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Figure  lb.  Zero-crossings 

Figure  1.     An  initial  function  ^i  and  the  zero  crossing  pixels  in  v,  t,  where  v,  j  blurs  the  ini- 
tial data  by  scaling  in  k. 
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Blurred  data 


Figure  2.  The  reconstructed  fui«;ti<^p  ^-  u^ij^  the  formula  from  Proposition  2.  The  data 
represents  the  result  of  blurring  the  g,  data  rcstrfcted  to  the  central  positive  interval,  zero  ex- 
tended elsewhere,  to  a  level  k  above  the4op  of  the  zero  crossing  for  v,  ^^  shown  in  Figure  1 
above.  The  reconstruction  uses  only  information  about  v,jt  along  the  zero-crossing,  and  is 
nearly  exact  to  machine  precision. 
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g 
Reconstructed 


Figure  3.  An  attempt  at  T's-ranjtru^ti.ngjj^/^^i^g^e  data  />,  from  Figure  2.  The  attempt  fails 
because  the  matrix  equation  Kg  =  0  relating  the  g  data  to  the  b  data  is  poorly  conditioned, 
even  though  there  are  mauy  more  b  values  tlian  g  values.  Thus  the  small  errors  due  to 
round-off  in  representing  the  i>,  data,  are  ma^ti'icd  when  reconstructing  g,.  The  singular 
value  decomposition  software  in  Crcvc  Meier's  "ni<itlab"  package  was  used. 


4.   Comments 

Zero-crossings  of  scale-space  filtered  data  seems  like  an  unlikely  form  of 
representation  of  data.  The  results  presented  here  suggest  that  even  when  supple- 
mented with  gradient  data  along  the  zero-crossings,  the  representation  is  still 
unlikely.  However,  the  instability  of  the  representation  does  not  completely  deny 
its  utility,  since  it  might  happen  that  the  classes  of  functions  mapping  to  similar 
representations  share  properties  essential  for  interpretation.  A  required  step  in  the 
validation  of  the  utility  of  a  representation  is  an  analysis  of  the  invariant  properties 
of  signals  that  yield  similar  representations.  A  necessary  condition  is  that  attempted 
reconstructions  differ  from  originals  in  unessential  ways,  (from  the  standpoint  of 
interpretation).  The  methods  outlined  here  should  prove  useful  in  verifying  or 
disproving  this  necessary  aspect  of  establishing  a  viable  representation. 
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